Peierls Instabilities in Quasi-One-Dimensional Quantum Double- Well Chains 
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Peierls-type instabilities in quarter-filled (n — 1/2) and half-filled {n — 1) quantum double-well 
hydrogen-bonded chain are investigated analytically in the framework of two-stage orientational- 
tunnelling model with additional inclusion of the interactions of protons with two different optical 
phonon branches. It is shown that when the energy of proton-phonon coupling becomes large, the 
system undergoes a transition to a various types of insulator states. The infiuence of two different 
transport amplitudes on ground states properties is studied. The results are compared with the 
pressure effect experimental investigations in superprotonic systems and hydrogen halides at low 
temperatures. 
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I. INTRODUCTION 

It is long been known that the formation of hydro- 
gen bonds between molecules or ionic groups is responsi- 
ble for a drastic changes in a wide variety of entire sys- 
tem properties such as structural x)hase transformations 
and proton ordering phenomenaO'D. In addition, proton 
transport phenomena in H-bonded materials and superi- 
onic properties discovered in some hydrogen-bonded crys- 
tals (for example, M3lI(A04)2 class where M=Rb, Cs, 
NII4; A=Se, S) are related closely to hydrogen-bonded 
network rearrangement. On heating these crystals trans- 
form into superionic conducting phase with statistically 
disordered hydrogen-bonded network (Fig. 1(a)). The 
protons can migrate through the two-dimensional con- 
ducting planes with low activation energy (~ 0.1 eV). In 
this case protonic conductivity increases significantly ti 
the value about 0.1 n~^- cm~^. It is generally accepted! 
that the two-stage conduction mechanism is required to 
sustain proton transport. The intrabond proton tun- 
nelling along the hydrogen bridge is connected with the 
transfer of ionic positive and negative charged defects, 
whereas the intermolecular proton transfer due to reori- 
entations of molecular group with proton leads to the 
breaking of the hydrogen bond and creation of a new one 
between another pair of molecular complexes. It should 
be noted that the formation of the hydrogen bridge in- 
duces the distortion of groups involved in hydrogen bond 
towards the proton that results in the shortening of the 
bonds. By this means the protonic polaron is localized 
between distorted ionic groups in the low-temperature 
ferroelastic phases, giving rise in this case to the dimer- 
ized structure. As has been shown in Ref. g, the small- 
radii polaron is formed due to the strong coupling of pro- 
ton with optical stretching vibration modes of the oxygen 
ions. It is evident that such an transformations from the 
superionic phase occurring in systems on cooling have the 
mixed (displacive and order-disorder) character. Theo- 
retical investigations of various ferroelectric-type order- 



ings in hydrogen-bonded systems have been based gener- 
ally on pseudospin Ising-type models with additional in- 
cluding of the pseudospin-phonon interactions to describe 
the coupling of protons with lattice vibration modes. In 
particular, the quantum double-well chain with quar- 
tic symmetric double-well potential has been used to 
model the transition from the symmetry-broken to the 
symmetry-restored ground state in hydrogen halides HX 
(X=F, Br, C1)eI which consist of hydrogen-bonded chains 
with weak interchain coupling. The dynamics of both 
ionic and orientational defects created by the rotations 
of molecular groups in hydrogen halides has been stud- 
ied inlhe framework of classic approach based on soliton 
modcH. 

It must be emphasized that taking into account the 
two-stage transport mechanism renders the pseudospin 
formalism unsuitable for the proton subsystem descrip- 
tion since the number of protons can differ from the num- 
ber of possible (virtual) hydrogen bonds and the proton 
occupancy of each bond in principle can be other than 
unity due to reorientational hopping and consequent fea- 
sibility of proton migration along the chain. Such a situ- 
ation, as an example, is observed in superionic materials 
of M3H(A04)2-type which transform oi|LCOoling into di- 
electric state with dimerized structureEl. It should be 
noted that such type of transitions to dielectric states 
is reminiscent of that of electronic systems in which the 
Peierls instabilities are observed. There have been many 
works to study the metal-insulator Peierls transitions in 
electron-phonon systems whiclp-|are unstable against the 
electron-phonon interactionaa'ta. It is common knowl- 
edge that the Peierls instabilities occur with the forma- 
tion of Peierls gap at fc = zLkp {kp is the Fermi level)in 
the electronic energy band that connected with the elec- 
tronic charge density waves condensation and structural 
lattice distortion modulations with q = 2kp. The ap- 
pearance of insulator state together with the structural 
transformation can be modelled in the framework of the 
Holstein electron-phonon model without additional in- 



eluding the anharmonic terms in the lattice potential. 

Recent investigations of Peierls transitions in electron- 
phonon systems have prompted us to study similar effects 
in several hydrogen-bonded solids. On the one hand, par- 
allel sequences of plains (001) are formed in supcrionic 
state of M3H(A04)2 crystals. These hexagonal conduct- 
ing plains consist of AO4 groups connected by virtual 
hydrogen bonds (see Fig. 1(a)). In the low-temperature 
phases the frozen-in hydrogen bonds with only one index 
/ (/ = 1,2,3) form well-defined sequences of dimers that 
involves the appearance of the parallel dimcrizcd chain 
arrays consisting of ionic groups linked by the /th hydro- 
gen bond (see Fig. 1(b)). To analyze the influence of the 
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FIG. 1. 
M3H(A04) 



(a) Hydrogen-bonded network in (001) plane of 
)2 crystal group; the solid lines indicate the possible 
type of dimerized structure which can appear with (/ — 3)th 
H-bonds frozen-in. (b) Structure of H-bonded dimer in one 
of low-temperature dimerized phases. 

proton-ionic group displacements coupling we consider 
the simplified model, namely the quasi-one-dimensional 
quantum double-well chain along one of the proton path- 
ways (for instance, the virtual hydrogen bond sequence 
. . . — 1 — 3 — 1 — 3 — . . .). As an initial step, we ne- 
glect of the interproton repulsion effect that is justified 
for the low proton concentration (in our case each pro- 
ton is averaged over every three virtual hydrogen bonds) . 
However, we take into account the possibility of proton 
exchange between our selected chain and surrounding. 
On the other hand, besides these superionic compounds 
we also analyze in this work the influence of ionic group 
displacements on the proton subsystem behavior in quasi- 
one-dimensional solid hydrogen halides. We reveal possi- 
ble symmetry-broken phases with proton charge dispro- 
portionalities coming from a Holstein coupling to AO4 
ionic groups or X atoms. We compare our conclusions 
with the results oL-the pressure effect theoretical stud- 
ies in M3H(A04)2EJ and hydrogen hahdcsB'E3. Although 
the first step of our analysis consists of the quasi-one- 
dimensional chain study, we believe our results can also 



be relevant for other hydrogen-bonded materials. 



II. DESCRIPTION OF THE MODEL 

The object of our consideration is the chain shown in 
Fig. 2(a). However, to avoid the geometric complexities 
introduced by the kinks in such an zig-zag chain, we con- 
sider in our model linear chain (see Fig. 2(b) where two 
neighboring chains are shown). The process of the pro- 
ton transfer in the double-well H-bond potential is rep- 
resented as the quantum tunnelling between two proton 
states with intrabond transfer integral f2o 
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where c^, cii, denote proton creation and annihilation 
operators in the position (/, v = a,b) of the chain. Be- 
sides that, we describe the interbond reorientational pro- 
ton hopping in two-level approximation as the quantum 
tunnelling effect with hopping amplitude fi^ 
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In this way in the framework of oricntational-tunnelling 
model proposed in Ref. |T^ the two-stage proton migration 
mechanism can be considered as the sequential migration 
of the ionic and orientational defects. 

As far as such a double-well chain is just a structural 
component of the system we also admit a possibility of 
proton exchange between the chain and surroundings by 
considering the system thermodynamics in the frame- 
work of the grand canonical ensemble with inclusion of 
the proton chemical potential 
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which is to be determined at the given proton concen- 
tration in the chain from corresponding equation for the 
chemical potential. 

Our main interest is to analyze the infiuence of the 
longitudinal optical ionic group vibration modes on the 
proton subsystem ground state. However, it was noted 
in Ref. n4 that the interactions between protons of the 
neighboring chains can lead to appearance of three- 
dimensional ordering. The more detailed analysis of the 
interchain proton interaction effect in this model together 
with the determination of stability conditions for the ex- 
istence of the phases with different orderi|nff type at finite 
temperatures will be presented elsewherdl3. We consider 
the anti-phase stretching vibration mode which causes a 
change of H-bond length in chain as indicated in Fig. 2(b) 
by solid arrows. Besides that, we also take into account 
the optical in-phase vibrations of ionic groups in chain 
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FIG. 2. (a) Zig-zag hydrogen-bonded chain in hydrogen 
hahdes, arrows indicate the possible path of proton migration 
along the chain, (b) Simplified model chains, the anti-phase 
and in-phase displacements of ionic groups identified by solid 
and dashed arrows. 



forms two energy bands with the bandwidth Ae = 
fio + ^fi — 1^0 — ^b\- The energy gap in this case is 
Aafc — 2|ilo — ^i?|- Eliminating one of the elementary 
transport process by setting hopping amplitude f^o = 
or rj/j = we can see that both the energy bands degen- 
erate into two energy levels and the quantum fluctuations 
between these two system states could be derived. It is 

clear that in the case when ?i = ;^ X) ^iv ~ 1 (one pro- 
lix 
ton is averaged within the bond) the lower band is filled 
and the chemical potential n is centered between bands - 
thus the material is in dielectric state. Such an situation 
can be observed in hydrogen halides. However, for n ~ h 
only half of the lower band is filled and this corresponds 
to the case of protonic conductor that occurs for example 
in superionic phases of superprotonic crystals. 

We will discuss afterwards the consequences of the 
proton-phonon coupling effect focusing on the analysis 
of the two physically different cases n — ^ (1/4-filled 



two-band model) and n 
model) . 



1 (half-filling case in two-band 



which induce their displacements with respect to sur- 
rounding chains as identified in Fig. 2(b) by dashed ar- 
rows. The coupling to the first type of displacements 
leads to the equal change of the potential well {I, a) and 
{I, b) depth within the H-bond 
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whereas the coupling of protons to another optical mode 
causes the difference of these potential minima depth 



'^Tl^ {q){nia ~ ni^i^b)(bq,2 + &^ 
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Here r/^^(g) = -2igiyfn/2MNuji{q) s\n\qdeyiY>[iq{l + 

l/2)d] and t/^^((j) = g2^/fij2MNuJ^ expWd] where 
gi and g2 are corresponding coupling constants, M is 
the effective ionic group mass, N denotes the number of 
hydrogen bonds in chain and d is a lattice spacing. Fur- 
thermore we take dispersionless approximation for the 
phonon frequencies: uii{q) = Wi and ^2(9) — ^2 and 
assume the harmonic approximation for the lattice vi- 
bration energies 
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First of all let us consider the case of the isolated chain 
without coupling to the phonon bath. Since the Hamil- 
tonian (0)-(m) can be exactly diagonalized, the proton 
energy spectrum 
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III. BROKEN-SYMMETRY SOLUTIONS 

A. Case n = ^ 

Let us now focus on the case of quarter filling when the 
half of the lower proton band is filled (one proton per two 
bonds). Then the macroscopic condensed phonau^tate 
is predominantly stabilized at q* = 2kF — 7r/dQ'Ej and 
is characterized by the expectation values of the phonon 
creation and annihilation operators 
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where A denotes the introduced distortion order param- 
eter which should be determined from the stationarity 
conditions of the free energy. Since the condensation of 
displacements (g) leads to the unit cell doubling, using 
the Fourier transformation Ci^ii\ - 



where the index i = {+,— } denotes {I — 2m) or 
{I = 2m + l)th cell, the Hamiltonian in condensed state 
with static periodic distortions (@) (adiabatic treatment) 
is given by 
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where Eq =r-Lhgi)'^ /2AI{'hu!i)'^ is well known from po- 
laron theoryllj protonic polaron binding energy which 



appears in the expression for the elastic energy per H- 
bond i^ and A = M^.JhpMZj[ = AA^/Eohuji/gi. 
The similar result can be obtained when we consider the 
second type of the ionic group displacements, in this case 
{Bq,2) = {bq,2 + blLq2) = —^NSq^q* and the Hamiltoniau 
in condensed state has the form similar to (P) with A — > 
A' = 2y/h/2Muj2 and Eq -^ E^ = {hg2f / 2M {huj2f . 
Since the inclusion of the coupling to the second phonon 
mode leads merely to renormalization of the binding en- 
ergy Eq in the Hamiltonian, further we focus on the anal- 
ysis of (|9|) with only one type of displacements taken into 
account. Introducing the double-time one-fermion diag- 
onal Green functions one can get rigorously the density 
of proton states 
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i?2 = V4^o(A2 + f^?j) - (tl - £2)2 

and the following notations arc introduced: ti = Vlf^ 
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A^ and 8(x) is the Heaviside 



step function. The expression for the ground state energy 
can be obtained easily from (0) and (UQ): 
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n%cos^ kd. (12) 



To determine the stable phase the equation dF/dA = 
should be solved. It appears that this equations has 
besides A = 0, a nonzero additional solution A 7^ 
for gi > gp where gp is the crossover proton-phonon 
coupling strength. The solution A 7^ corresponds 
to the global minimum of F and, as a result, implies 
the structural distortion stabilization with the amplitude 
ui = y/h/2MNuji{Bq>) = ^(-1)' (see Fig. 3(a)). Let 
us discuss the proton position average occupancies on the 
bonds and the band structure. At gi — gp each proton 
band splits into two subbands 



-'a{+/-)( 



{k)^T\ti+2noJA'2 + nlcos^kd, (13) 



£b{+/-){k) =±yh- 2f7oyA2 + rj2^cos2 kd 

as shown in Fig. 3(b) where the proton density of states in 
the disordered and dimerized phases is represented. The 



Peierls energy gap between either of the two (lower and 
upper) subbands Ai — \Jt\ -I- 2rio^ ^ v ii ~ 2rio^ ~ 
2fioA/-\/^o + ^R ^"^^ tends to zero for f^o -^ 0. In 
this case A 



\p and tends to zero for 



rijj and the phase transition 
(change in the nature of the ground state) occurs when 
the localization energy Eq ^ {gXY — \^r- The energy 
gap between the second and third subbands increases at 
.91 > 9P 



Aab^2\lti^2^a\IA^ + n\. 



The proton chemical potential /i is centered between 
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FIG. 3. (a) Distortion parameter A' = j-j- as a function of 

proton-phonon coupling g\ for f^o = 0.5; inset: dependence of 
average proton accupancies on g\ for CIh — 0.14. (b) Proton 
density of states p(e) = j^ (e = j-^), dashed and dotted 
curves indicate the cases of A' — 0.5 and A' — 0.0 respec- 
tively. 

two lowest subbands with further increasing of gi > gp 
that points to the insulator state appearance. We see 
from inset in Fig. 3(a) that the distortion stability is 
accompanied by the formation of the proton charge- 
density-wave state in which (ti/q) = {nn,) = j{l + (—1)') 
that means the forming of dimerized structure as shown 
in Fig. 4. Consider further the ground state phase di- 
agrams (51 = gi/huji, 0,0 = flo/huji) and (51, Op = 



^^^9-^ 



FIG. 4. Dimerized structure which appears in the case of 
quarter-filled chain. 



ilfj/huji) represented in Fig. 5. We see the strong influ- 
ence of the amplitude Qr on the dimerized state stabil- 
ity. The increasing of ^Ir suppresses dimerization. At 
riji — > (without reorientational hopping) the system 
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creases and approaches the critical value g* = ^/flIi/2, a 
transition from the dimerized to the uniform state occurs 
with pressure. This effect appears due to the more weak 
proton-phonon coupling and, as a result, to the tendency 
of the proton delocalization in chain. 



B. Case n = 1 

Let us discuss another case when one proton in ajuerage 
is placed in the bond. According to Peierls theoryO such 
a system is very susceptible towards lattice modulation 
at g* = 0. It should be noted that in this case only 
the second type of optical vibrations (interchain mode) 
contributes to the lattice distortions condensation. The 
Hamiltonian in condensed phase has the form 



H = 



y^ I (-A* + A)nfea - (a* + A)nfeb 

k 
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FIG. 5. Ground-state phase diagrams (a) (gi, Qq) and (b) 
(gi, Qr). The notations PD and PU denote the dimerized 
and uniform phases respectively. 

is brought immediately into the dimerized state. Only 
for the finite values of Qr the uniform disordered phase 
begins to appear and the "metal" -insulator transition oc- 
curs. 

It is necessary to mention that the hopping amplitudes 
ilo ^-nd Qr depend strongly on external pressure. In 
particular, the Jig value decreases with pressure-that de- 
duces from quantum mechanical calcuLartionaill as well 
as from the experimental measurementsEJ. This is asso- 
ciated with the shortening of the distance between two 
potential minima {I, v) in the bond. Thus we can make 
a conclusion about pressure effect on the system state 
from the diagrams shown in Fig. 5. Using the obtained 
ind values for parameters VIr, gi and wi (ytR/huji sa 0.14 
and h?'g\/2M(Tiu}if « 3.8) we reveal that the dimer- 
ized state is always stable at T=0 under pressure for 
this set of parameters. It is interesting that the simi- 
lar picture has been observed in M3H(A04)2 materials 
from experirneatally measured baric dependencies at low 
temperatureali. Nevertheless, we notice that as gi de- 



In this case the density of proton states 
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points to the two-band structure 



e,(fc) = ±VA2 + |ifeP 
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with the Peierls energy gap Aab = 2-wA2 -I- (fio — ^rY- 
The chemical potential always is centered between two 
bands, i.e. fi = 0. We present the equation for deter- 
mination of A 7^ which follows from the stationary 
condition of F: 
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The nonzero solution, which appears for g2 > gp, corre- 
sponds to the formation of a proton charge density wave 
in chain together with the distortions ui = A/252 stabi- 
lization (see Fig. 6) . The typical dependencies of the pro- 



FIG. 6. Broken-symmetry structure which appears in the 
case of half-filled chain. 

ton position average occupancies {niu) are represented in 



Fig. 7. It is interesting that the system now is invariable 
with respect to the interchanging f2o '^ ^R- Thus it is 
sufhciently to analyze the system behavior as a function 
of CIq for instance, with the given fixed value of ^Ir. The 
ground state phase diagram (32 — 32/^1^2, CIq — flQ/huj2) 
(see Fig. 8) differs essentially from the case n — ^. The 
phase equilibrium curve has the specific salient point at 
r^o = ^R- The drastically decrease of gp in the vicinity 
of r^o = ri/j is connected with the fact that the transfer 
anisotropy |r2o~^i?| forms the additional transverse field 
which competes with the ordering stabilization process. 
This anisotropy is vanished at fig = ^R. that leads to the 
lowering of the crossover proton-phonon coupling energy 




FIG. 7. Average proton occupancies as a functions of 32 
for Qq — 0.8; bold and thin curves indicate the cases when 
CIr = 0.5 and CIr = 2.5. 

gp required for the ordering stabilization. The interpre- 
tation of the diagram (^2, ^0) with respect to the pres- 
sure effect is very interesting. The second-order transi- 




12 3 4 5 

FIG. 8. Ground-state phase diagrams (g2, S7o). The no- 
tation PO indicates the symmetry-broken phase with proton 
ordering on the hydrogen bonds; inset: the region flo ~ Qr 
in more details. 

tion from the uniform to ordered state occurs under pres- 
sure (with fio decrease) for 32 > 52 — V^hi/2. However, 



in the region (72 < 52 ^^^ additional reentrant transition 
from the symmetry-broken to the uniform state appears 
(see Fig. 9). In this case the region of symmetry- broken 
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FIG. 9. Distortion parameter as a function of D,r for 
§2 = 0.32 and flo = 0.5; inset: corresponding dependencies of 
average proton accupancies. 



ase. We notice 
'^^ and the results 



phase equilibrium narrows with g- 
that the first-principle calculationsE3 
of Monte Carlo simulationaj in quasi-one dimensional 
hydrogen halides show a transition from the symmetry- 
broken phase shown in Fig. 6 to the uniform symmetric 
phase under pressure at the low temperatures. Thus our 
results in the vicinity of fla w Up and for considerably 
weak proton-phonon coupling 52 < 32 ^^^ ^^ qualitative 
agreement with the conclusions of Refs. HJ13 confirm- 
ing a proper treatment of the quantum effects in these 
hydrogen-bonded materials. 



IV. CONCLUSIONS 

In the present work the lattice effect on the ground 
state properties of the quantum quasi-one-dimensional 
hydrogen-bonded chain is analyzed in the framework of 
the two-stage orientational-tunneling model. The inter- 
action of protons with two different types of surrounding 
ionic group optical displacements is considered. We show 
that when the proton-phonon coupling energy becomes 
large, the system undergoes a transition from disordered 
to broken-symmetry phases. The different cases of proton 
concentration have been analyzed: n = 1/2 and n = 1. 
It is shown that in the first case the Peierls transition 
to the dimerized phase occurs, whereas in the second 
one we obtain a transformation into the proton-ordered 
state. The infiuence of two different transport amplitudes 
on ground states properties is also studied. We compare 
our ground-state phase diagrams with the pressure effect 
experimental investigations in superprotonic systems and 
hydrogen halides at low temperatures. 
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